Effects of isoelectronic Ru substitution at the Fe site on the energy gaps of optimally 

F-doped SmFeAsO 



D. Daghcro, 1 M. Tortcllo, 1 G.A. Ummarino, 1 V.A. Stcpanov, 2 
F. Bcrnardini, 3 M. Tropeano, 4 M. Putti, 5 ' 4 and R.S. Gonnelli 1 

1 Dipartimento di Scienza Applicata e Tecnologia, Politecnico di Torino, 10129 Torino, Italy 
2 P.N. Lebedev Physical Institute, Russian Academy of Sciences, 119991 Moscow, Russia 

'Dipartimento di Fisica, Universitd di Cagliari, 0904-2 Monserrato ( CA ), Italy 
* Dipartimento di Fisica, Universitd di Genova, via Dodecaneso 33, 16146 Genova, Italy 
5 CNR-SPIN, Corso Perrone 24, 16152 Genova, Italy 

We studied the effects of isoelectronic Ru substitution at the Fe site on the energy gaps of optimally 
F-doped SmFeAsO by means of point-contact Andreev reflection spectroscopy. The results show 
that the SmFei_ x Ru x AsOo.85Fo.is system keeps a multigap character at least up to x — 0.50, 
and that the gap amplitudes Ai and A2 scale almost linearly with the local critical temperature 
Tg. The gap ratios 2Ai/kBT c remain approximately constant only as long as T c > 30K, but 
increase dramatically when T c decreases further. This trend seems to be common to many Fe- 
based superconductors, irrespective of their family. Based on first-principle calculations of the 
bandstructure and of the density of states projected on the different bands, we show that this 
trend, as well as the T c dependence of the gaps and the reduction of T c upon Ru doping, can be 
explained within an effective three-band Eliashberg model as being due to a suppression of the 
superfluid density at finite temperature that, in turns, modifies the temperature dependence of the 
characteristic spin-fluctuation energy. 

PACS numbers: 74.50. +r , 74.70.Dd, 74.45. +c 



INTRODUCTION 

The discovery of Fe-based superconductors (FcBS) [l[ 
with T c as high as 55 K has shown that cuprates no longer 
represent the only class of high-T c compounds. One of 
the reasons of the great excitement in the scientific com- 
munity and of the impressive amount of work produced 
up to now is certainly that these materials give the op- 
portunity to study high-T c superconductivity in different 
systems, in the hope to enucleate its key elements. The 
parent stoichiometric compounds of most FcBS are not 
superconducting (with few exceptions, like LiFcAs and 
LaFePO) but display a metallic behaviour (as opposed 
to the Mott insulating state of cuprates) and feature a 
long-range antiferromagnetic (AFM) spin-density-wave 
(SDW) order. Superconductivity appears upon doping 
and, in some systems, also by applying pressure; how- 
ever, the order of the transition between magnetic and 
superconducting phases seems not to be universal though 
increasing evidences are being collected of a region of co- 
existence of superconductivity and magnetism. Contrary 
to cuprates (where the superconducting region in the 
phase diagram is dome-shaped and the maximum T c cor- 
responds to a well-defined "optimal" doping) in FcBS su- 
perconductivity sometimes appears with T c already very 
close to the maximum and shows a weak doping depen- 
dence in a broad doping range. A central feature of FeBS 
- which is tightly connected to the origin of superconduc- 
tivity according to the most widely accepted theories - 
is their multiband character. They feature indeed two or 
three hole pockets around the V point of the first Bril- 



louin zone and two electronpockets at the M point (in 
the folded Brillouin zone) @, Q. In 1111 compounds, 
all the relevant Fermi surface sheets are weakly warped 
cylinders parallel to the k z axis (as expected in a layered 
material) while a greater degree of three-dimensionality 
is observed in 122 compounds. These multiple bands and 
their almost perfect nesting in the parent compound ex- 
plain the AFM instability The weakening of the nesting 
induced by doping instead leads to spin fluctuations that 
would act as the glue for the formation of Cooper pairs. A 
spin-fluctuation mediated pairing would be mainly inter- 
band and would favour the opening of superconducting 
energy gaps of different sign on different Fermi surface 
sheets, the so-called s± symmetry 0]. Though many 
theoretical and experimental results support this theory 
0, Q there is not, up to now, a definitive proof of such 
a picture. Things arc further considerably complicated 
by the fact that the electronic banstructure is very sensi- 
tive to some fine structural parameters, like the Fe-As-Fe 
bond angle and more particularly the height of the pnic- 
togen atom (IiAs) above the Fe layer. Possibly because 
of this sensitivity, in many situations the gap structure 
of FeBS can vary considerably within the same system, 
giving rise to line nodes, point nodes, deep gap minima 
etc. [7H10I| . In 1111 compounds, Ila s has been proposed 
as a switch between high-T c nodeless superconductivity 
and low-T c nodal superconductivity 
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In the effort to discriminate the effects of different pa- 
rameters on the superconducting and magnetic phases 
of FeBS many different chemical substitutions have been 
performed. The main effect of aliovalent substitutions is 
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to dope t he p arent compound with charge, either elec- 



trons [12J, |13[ or holes [14j thus allowin g to explore the 
phase diagram. Isovalent substitutions IB- HI, instead 
were tried to modify the lattice structure, create "chem- 
ical pressure" , introduce disorder (acting as magnetic or 
non-magnetic impurities) etc. A further degree of free- 
dom is the site of substitution, that can reside either in 
the spacing layer [l4[ or in the active one containing Fe 
IB" HI, which is possibly directly involved in the mag- 
netic pairing via spin fluctuations. 

Here we report on point-contact Andreev-rcflcction 
spectroscopy (PCARS) measurements performed in 
SmFei_a;Ru x AsOo.85Fo.i5 with x ranging from to 0.50. 
The considerable decrease of T c in this series of sam ples 
has been attributed to disorder in the Fe sub-lattice [15( 
and/or occurrence of a short-range static magnetic or- 
der [l9j]. PCARS results clearly indicate the presence 
of a multigap character at all the investigated levels of 
Ru substitution. The superconducting gaps decrease ap- 
proximately linearly with the local critical temperature 
of the contact, but, even when the latter is reduced by 
a factor 5 with respect to the optimal value, they show 
no sign of nodes, either intrinsic or "accidental". For 
both gaps, the 2A/fcsT c ratio is rather constant down to 
T C A > 30if but then increases consistently below this crit- 
ical temperature. Comparison with other results in liter- 
ature indicates that many different FeBS fit in this trend, 
which suggests the possibility to study some properties 
common to different compounds in a single samples series 
that allows spanning a very wide range of critical temper- 
atures. Thanks to ab-initio electronic structure calcula- 
tions, the trend of the gaps as a function of has been 
reproduced within a minimal three-band, s± Eliashbcrg 
model. This model also takes into account the so-called 
"feedback" effect, i.e. the effect of the condensate on the 
antiferromagnetic spin fluctuations possibly responsible 
for the superconductivity in these compounds. The evo- 
lution as a function of T c of the temperature dependence 
of the condensate necessary to reproduce the experimen- 
tal data looks rather similar to that obtained from Lon- 
don penetration-depth measurements performed in other 
FeBS, particularly in the region of coexistence of super- 
conductivity and magnetism. This fact suggests that, in 
agreement with ref. (l9j . proximity of superconductiv- 
ity and magnetism in these samples might be one of the 
main reasons for the decrease of T c and for the observed 
behavior of the energy gaps. 



EXPERIMENTAL DETAILS 

The polycrystalline SmFei_ x Ru x AsOo.85Fo.i5 samples 
were synthesized as described in Ref. [15j. The starting 
mixture of fine powder of SmAs and 99.9% pure Fe203, 
Ru02, FeF2, Fe and Ru was pressed in pellets and then 
put through a two-step reaction process involving a first 
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FIG. 1: Lattice constants for SmFei-^Ru^AsOo.ssFo.is at 
different Ru content x. Lines are only guides to the eye. 
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TABLE I: Resistive critical temperatures and residual resis- 
tivities (defined as in ref. 0) for SmFei 

— X RUajAsO0.85F0.15 

samples at different Ru contents. T c and po for the samples 
with x — 0.25 and x = 0.50 are different from those reported 
in ref. [H. 



heating to 450° and a second heating to 1000-1075°. X- 
ray diffraction analysis showed small amounts of SmOF 
(up to 6 %) in the final samples. Figure [T] shows the de- 
pendence of the lattice constants a and c on the Ru con- 
tent x, indicating that Ru substitution for Fe is effective. 
Resistive critical temperatures and residual resistivities 
for the samples used in this work are reported in table [I] 
The samples with x = 0.25 and x — 0.50 have higher T c 
and much improved transport properties (namely, resis- 
tivity, magnetoresistance and Hall mobility) than those 
reported in ref. |15j | for the same doping contents, even 
though they were prepared in the same way. The possi- 
ble reason of this difference is under investigation. In any 
case, these samples were particularly suited for PCARS 
measurements, since the longer mean free path makes it 
easier to attain the spectroscopic conditions, as explained 
below. 

Point contact spectroscopy is a local, surface-sensitive 
technique and it is therefore necessary to avoid any sur- 
face degradation or contamination. The samples were 
thus always kept in dry atmosphere, and broken to ex- 
pose a clean surface prior to point-contact fabrication. 
The point contacts were made by putting a small drop of 
Ag paste on that surface, as described elsewhere [2(J 21 
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With respect to the standard "needle-anvil" technique, 
this configuration ensures a greater mechanical and ther- 
mal stability of the contacts and also allows the whole 
mounting for point contact to be hermetically closed in 
the cold head of the cryogenic insert thus avoiding any 
exposition to air and moisture during the transfer from 
the glove box (where the point contacts are fabricated) 
to the cryogenic environment. Although the Ag drop has 
a diameter of at least 50 fim, the real electric contact oc- 
curs only between some of the Ag grains and the sample 
surface. The true contact is thus the parallel of sev- 
eral nanoscopic junctions that can well be in the ballistic 
regime (i.e. have a radius smaller than the electron mean 
free path). In Ref. [l5[ a rough evaluation of the mean 
free path in SmFei_ x Ru x AsOo.85Fo.i5 gave t = 3 — 10 
nm without any clear dependence on the Ru content. In 
the cleaner samples with x — 0.25 and x = 0.50, the 
same evaluation gives I ~ 7 nm and I ~ 20 nm, respec- 
tively. Such small values of the mean free path make 
the fulfillment of the ballistic condition a <C £ (where a 
is the contact radius) be very difficult to achieve. For 
instance, with these values of £ and the residual resistiv- 
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ities taken from table HI the Sharvin equation 22| would 
require resistances of the order of several kfi for the con- 
tact to be ballistic. The typical experimental resistance 
of the contacts is instead in the range 10 — 100f2. Indeed, 
many of the contacts were not spectroscopic or showed 
heating effects. A large number of measurements was 
then necessary to achieve a relatively small number of 
successful measurements. All the results reported here, 
except those shown in figlH are thus referred to the small 
fraction of contacts that do not show heating effects and 
gave a clear Andrccv-rcflection signal. In these cases, 
the existence of many parallel nano junctions can be in- 
voked to reconcile the actual contact resistance with the 
requirement of ballistic transport [2l|. In some cases, 
the Sharvin condition was fulfilled at low temperature 
but broke down on increasing the temperature because 
of the decrease in the mean free path. In these cases, 
the values of the gaps at low temperature can be taken 
as meaningful anyway, though their temperature depen- 
dence and eventually the value of the local critical tem- 
perature can be slightly affected by the non-ideality of 
the contact. 



RESULTS AND DISCUSSION 
Point-contact Andreev-reflection results 

Figure reports the temperature dependence of the 
raw conductance curves (obtained by numerical differen- 
tiation of the I — V characteristics) of one of the contacts 
that did not show any anomaly. The curves were mea- 
sured in the x = 0.10 sample, and the normal-state resis- 
tance of the contact is around 50 ft. The curves show the 
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FIG. 2: Temperature dependence of the raw conductance 
curve of a point contact on the x — 0.10 sample. The normal 
state resistance is Rm = 50J1. The inset shows the curves at 
4.2 K and at 21.0 K in an extended voltage range, to highlight 
the excess conductance persisting up to about 100 mV. 



typical features already observed in SmFcAs(Oi_ x F x ) 
23, in LaFeAs(Oi_ x F x ) [H and in other 1111 com- 
pounds. In particular, they feature clear maxima related 
to a presumably nodeless gap, shoulders suggestive of a 
second larger gap and additional structures that, as re- 
cently shown [2j4J, can be explained as being due to the 
strong electron-boson coupling. The excess conductance 
at high voltage, extending up to about 100 mV (see inset) 
is also typical of these systems [24| . The temperature at 
which the Andreev-reflection features disappear and the 
conductance becomes equal to the normal-state one is the 
local critical temperature of the contact, or Andreev crit- 
ical temperature T^. As shown in figf2]this temperature 
is easy to identify in spectroscopic contacts because it 
also marks the point where conductance curves recorded 
at slightly different temperature start to be superimposed 
to one another (here the curves at 21.0, 22.0 and 23.0 K 
coincide within the experimental noise). 

In contrast, Figure [3] reports two examples of conduc- 
tance curves that show, together with an Andreev sig- 
nal, deep and wide dips that, at low temperature, occur 
at energies comparable to those of the large gap. As 
shown elsewhere 21, 25[ these dips are likely to be due 
to the current becoming overcritical in the region of the 
contact and prevent a proper determination of the gap 
amplitudes. On increasing temperature, they move to- 
ward lower voltage (due to the decrease of the critical 
current) causing an apparent shrinkage of the Andreev 
signal, finally giving rise to a sharp cusp at zero bias. 

Going back to the case of ballistic contacts as in Fig[2~J 
the conductance at or just above can be used to nor- 
malize all the curves at T < T C A . In principle, a conduc- 
tance curve recorded at a given temperature should be 
normalized to the normal-state conductance at the same 



4 




-60 -40 -20 20 40 60 80 

Voltage (mV) 



40 -20 20 40 60 

Voltage (mV) 



FIG. 3: Temperature dependence of two raw conductance 
curves in case of non ballistic contacts in the x = 0.1 sample. 
Dip features characteristic of non-ideal conduction through 
the contact are clearly visible; the conductance minima are 
indicated by arrows. 
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temperature, but because of the very high upper critical 
field of these materials, the latter is not usually acces- 
sible, at least at low T. Using the normal state at 
to normalize all the curves is thus a somehow arbitrary 
choice but, as shown elsewhere [24j |. is anyway the one 
that preserves the weaker structures, i.e. those due to 
the large gap and, if present, those due to the strong 
electron-boson coupling. 

Figure 0] shows some examples of low-temperature, 
normalized conductance curves in samples with differ- 
ent Ru content. Some important points are immediately 
clear by looking at these curves. First, none of them 
displays zero-bias peaks, and the same happens in 100% 
of the spectroscopic contacts. This points towards the 
absence of line nodes even at the highest Ru contents, 
contrary to what has been observed in some other FcBS 
away from optimal doping 01- Second, all the 

curves show more or less marked double-gap features. 
Third, despite the very large range of doping, the width 
of the structures does not change very much (note that 
all the panels have the same horizontal scale). Thus we 
should not expect major variations in the gap values upon 
Ru doping. Fourth, the asymmetry of the normalized 
conductance curves for positive/negative bias - which is 
particularly strong in unsubstituted SmFeAs(Oo.85Fo.i5) 
[23( - seems to be reduced by Ru doping. As a matter 
of fact, it is clearly visible even at a first glance in the 
case x = 0.05, becomes discernible only while trying to 
fit the data in the cases x = 0.10 and x = 0.25 but al- 
most completely disappears for x = 0.30 and x = 0.50. 
The real origin of this asymmetry, which is common to 
most point-contact spectra in Fe-based superconductors, 
is not completely clear yet, th oug h it has been recently 
ascribed to the Seebeck effect [27|. Preliminary Seebeck 



FIG. 4: Some examples of experimental conductance curves, 
after normalization (symbols), of point contacts on samples 
with different Ru content. The curves were all measured at 
4.2 K. Solid blue (red) lines represent the best fit of the right 
(left) side of the experimental curves obtained within the two- 
band 2D generalized BTK model [H, [2(|. The asymmetry 
decreases on increasing the Ru vanishing completely at x — 
0.50. The fitting parameters are indicated in the labels. 



effect measurements performed in these samples show in- 
deed a considerable decrease of the Seebeck coefficient 
with increasing Ru content [28| . 

To extract quantitative information about the ampli- 
tude of the gaps from the conductance curves, they must 
be compared with suitable theoretical models. None of 
the models for single-band superconductivity can repro- 
duce the shape of the experimental curves of figHl In- 
stead, a two-band Blonder-Tinkham-Klapwijk model [29[ 
generalized to the 2D case [26[ and including a broaden- 
ing term (30l | is the minimal model that can be used in 
this case. For each band the parameters of the model are 
the energy gap A, the broadening parameter T and the 
barrier parameter Z. Then, being the total conductance 
the weighed sum of the single-band conductances, the 
last parameter is the weight w\ of band 1 (the weight of 
band 2 being consequently determined as 1 — wi).[45j It is 
true that Sm-1111 is not two-dimensional and thus a 3D 



model should be used; however, as shown elsewhere [24j, 
the latter is much more complicated and for any prac- 
tical purpose one can safely use the 2D one (especially 
when, as it is the case here, the gaps are nodeless). The 
lines in Fig |4] represent the best fit of the experimental 
data, and the labels indicate the relevant values of the 
gaps and of the weight of band 1 in the point-contact 
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conductance, w\ . To account for the residual asymmetry 
of the normalized curves, we chose to fit the positive- and 
negative-bias side separately (blue and red lines, respec- 
tively). As previously stated, for x > 0.30 the asymmetry 
is very small and the difference between the two fits can 
be no longer appreciated. 

Figure [5] shows two examples of how the normalized 
conductance curves evolve with temperature, and the re- 
sulting temperature dependence of the gaps extracted 
from the fit. The two cases shown refer to a lightly 
doped sample (x = 0.10 ) and to a heavily doped one 
(x = 0.50). The lowest-temperature curves show clear 
shoulders related to the larger gap, which become less 
and less discernible in the other curves (vertically off- 
set for clarity). In the x — 0.50 case, a dip structure 
is also seen to shift to lower energy on increasing the 
temperature, possibly giving rise to the small downward 
deviation of the temperature dependence of A 2 from a 
BCS-likc A(T) curve observed at high temperature. Al- 
though there is no reason to expect the gaps to follow 
a BCS-like curve, the effects of the dip do not allow us 
to discuss whether this deviation is intrinsic or is an ar- 
tifact due to the small mean free path of the samples. 
Incidentally, on the basis of recent calculations within a 
minimal three-band Eliashberg model j3l| one would in- 
stead expect the gaps to be greater than the BCS value 
in proximity of the critical temperature. 

Let us just recall here that the fitting procedure is 
generally not univocal, i.e. different sets of parameters 
can give almost equally good fits. Error bars in the insets 
to FigJS] indicate the spread of gap values resulting from 
different fits. 

Figure [6ja) shows the behavior of the gaps as a func- 
tion of the Ru doping x. The data are rather scattered 
but a general trend is anyway discernible. While the 
small gap Ai docs not vary sensibly on increasing the 
Ru content x, the large gap A 2 shows a rapid decrease 
from x = 0toa;=:0.10 and then remains approximately 
constant. This behavior is in rough qualitative agree- 
ment with that of the bulk T c reported in table fl] Since 
PCARS is a local probe, the scattering of gap values at 
the same composition is most probably due to slight in- 
homogeneities in the local doping content. As long as T c 
has a strong dependence on the doping content, i.e. up 
to x = 0.25, different point contacts on the same sample 
can thus provide different values of the gaps and of the 
local , i.e. the Andreev critical temperature. As a 
matter of fact if one plots the gaps as a function of 
as in figHfb), a roughly linear trend of both Ai and A2 
can be appreciated despite the fluctuations in their val- 
ues. It is worth reminding that the data reported here 
are already the results of a very careful selection aimed 
at eliminating all the questionable results, so that these 
fluctuations are not due to spurious effects that can be as- 
cribed to non-ballistic conduction, heating, or spreading 
resistance. As for the large gap, a large uncertainty was 
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FIG. 5: Temperature dependence of the normalized conduc- 
tance curves (symbols) of point contacts on the samples with 
x = 0.10 (a) and x — 0.50 (b) with the relevant two-band fit 
(lines). All the curves but the top ones are vertically shifted 
for clarity. The insets show the temperature dependence of 
the gaps Ai and A2 as extracted from the fit. The error 
bars indicate the spread of gap values obtained by different 
fits, when the other parameters (Ti, Zi and the weight wi are 
changed as well). 



also found in the starting compound SmFeAs(Oi_ x F x ) 
[23| (here represented by the vertical error bar on the 
point at the maximum T^) and was ascribed to the resid- 
ual degrees of freedom in the normalization process, to 
the asymmetry of the curves and to the fact that the 
features related to A2 are less sharp than those related 
to the small gap Ai. However, in the high-doping range 
(x = 0.3 and x = 0.5), depends very little on the 
Ru content and the asymmetry has almost completely 
disappeared; even large differences in local composition 
correspond to a small difference in T^. Therefore, the 
spread of gap values accompanied by a small spread in 
seems to indicate a lack of correlation between these 
quantities (as observed also in MgB2 with AI and Li co- 
doping [32J). 

Since the gaps show an overall linear trend as a func- 
tion of the local , it is particularly instructive to plot 
the gap ratios 2Ai/fcsT c as a function of T C A . This is 
done in figure[71 which also reports various other PCARS 
data in 1111 and 122 compounds. Only results show- 
ing nodeless order parameters are shown for consistency. 
It is clear that the Ru substitution in the optimally F- 
dopcd Sm-1111 allows spanning a wide range of criti- 
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FIG. 6: (a) Energy gaps as extracted from the two-band fit of 
the conductance curves, as a function of the Ru content x. (b) 
Energy gaps as a function of the local critical temperature 
of the various point contacts. Lines are gaps calculated within 
the three-band, s± Eliashberg theory (see text for details). 
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FIG. 7: Gap ratios 2Ai/fcsT c and 2A2/A:bT c as a function 
of T c as determined by point-contact Andreev-reflection spec- 
troscopy in Sm(Fei_ x Ru x )AsOo.85Fo.i5 (blue and red circles) 
and in various other 1111 and 122 materials. Solid lines are 
2Ai/fcsT c ratios calculated within the three-band, s± Eliash- 
berg theory (see text for details). Experimental data from 
literature are taken from ref. [M and references therein. 



cal temperatures, which not only covers but also extends 
the range of T c values of superconducting Fe-based com- 
pounds measured so far by PCARS. As already shown 
in ref. [24 |. the 2Ai/ksT c ratios start to increase be- 
low T c ~ 30 K. Surprisingly, PCARS data on this single 
sample series, namely SmFei-^Ru^ AsO0.85F0.15, feature 
basically the same behavior as those obtained from many 



other different nodeless FeBS of the 1111 and 122 fami- 
lies. 

These results appear to be in contrast to what reported 
in ref. [33j . where an opposite trend is suggested. How- 
ever, even in the aforementioned paper, several results 
reported for FeBS show 2A/fcsT c ratios which seem to 
increase with decreasing T c , particularly for the large gap. 
A definitive answer on the possibility of a universal trend 
of 2Ai/fcsT c vs. T c requires more experimental work, 
comparing results obtained with different techniques on 
samples of increasingly better quality. 



Electronic structure calculations 

In order to try to explain the observed PCARS data 
within the Eliashberg theory, we preliminarily performed 
electronic structure calculations. In particular, Sm- 
FeAsO electrons and holes density of states (DOS) have 
been obtained by ab initio calculations performed in 
the local density approximation to the density-functional 
theory (LDA-DFT) 



34| as implemented in the all 
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electron full-potential APW and local orbitals 
code Wien2k [37j |. APW has the advantage to treat ex- 
plicitly Sm 4f electrons within the valence band yield- 
ing state-of-the-art band structure dispersion quality. 
To simulate Ru substitution our calculations were per- 
formed in a tetragonal super-cell containing four formula 
units [Pma2 No. 28) where 25% Ru concentration was 
achieved by the substitution of one Ru out of four Fe, 
retaining the bulk symmetry in the defected cell. The 
conservation of symmetry allowed a reliable compari- 
son of doped and undoped band structures without fold- 
ing/unfolding mapping problems. Muffin tin radii of 2.3, 
1.9, 2.2 and 2.0 Bohr were used for Sm, O, Fe and As, 
respectively. Brillouin zone integration was performed 
with tetrahedrons on a 6 x 6 x 4 mesh [38| . Since we are 
interested in the superconducting paramagnetic phase of 
SmFeAsO, Fe spin polarization was not considered. The 
relevant band structure of this system is given by two 
electron and three hole bands as it is usual for 1111 iron- 
pnictides. Those are superimposed to the -£f Sm bands at 
the Fermi energy (Ep). To disentangle the contribution 
of the Fe 3d bands we used the so-called fat-band rep- 
resentation by projecting the wavefunction onto the Fe 
atomic orbitals and thus obtaining the band structures 
shown by dots in FiglU The size of dots is proportional to 
the d atomic character of the wavefunction on Fe atoms. 

Such a procedure allows to identify a number of 
parabolic Fe 3d bands around Ep. Since comparison 
with LaFeAsO is relevant, in Fig|S][a) we report the band 
structure obtained for the same crystal structure of the 
SmFeAsO compound where Sm has been replaced by La. 
In FigJSJa) we clearly see the usual set of five parabolic 
Fe d bands, two electrons (ei, 62) and three holes ones 
(hi, I12 and ha), all of them centered at the T point since 
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FIG. 8: Color online. SmRu^Fei-xAsO fat-band structure in the four formula unit super-cell. Dots are proportional to the Fe 
3d character of the wavefunctions. Parabolas show the results of a fit of the Fe 3d bands within the two-dimensional electron gas 
model. Upper and lower horizontal dashed lines are Ep with and without F doping, respectively, a) Reference band structure 
with La replacing Sm. b) Actual SmOFeAs band structure including the 4} Sm bands. The two light holes bands have similar 
dispersions and have been fitted with a single parabola hi. c) Same as in b) but in the case of SmFeo.7sRuo.25AsO. In b) and 
c) the parabolic fit of band 63 is shown only up to an energy close to the one where its hybridization with Sm 4f bands takes 
place. 



we are using a four formula unit super-cell. This find- 
ing is in agreement with previous LDA-DFT calculations 
provided wo consider that SmFeAsO structure was used 
here As for the band structure of the actual Sm- 
FeAsO, we find Fe d bands superimposed to Sm 4 f ones. 
This fact not only makes the interpretation of the band 
structure less easy but introduces an hybridization effect 
between Sm and Fe states. LDA bands in Fig|8](b) show 
that Fe electron bands in SmFeAsO are no more simply 
parabolic but hybridization introduces warping out of the 
central part of the Brillouin zone. Hybridization splits 
the electron band labeled e2 in FigJSJa) in two pieces 
named e 3 and ei in Fig|8](b). In undopcd SmFeAsO e3 
and ei do not cross Ep being the former too low in energy 
and hybridized with Sm 4f bands (see details of FigJS^b)) 
and the upper too high. Therefore, undoped SmFeAsO 
will have only one Fe-derived electron band e2 and three 
hole ones, the doubly-degenerate light-hole hi and the 
heavy- hole I12. The Fermi surface will be made of four 
nearly cylindrical sheets, one less than LaFeAsO. Such 
a finding seems to be in agreement with ARPES mea- 
surements on SmFeAsO where only one electron band is 
suggested j3§] . Since in the Eliashberg approach we are 
interested in the Fe d states DOS, we need to disentangle 
that contribution out of the total one that includes the 
Sm 4f bands. Our choice is to model the band structure 
as a superposition of five parabolic bands (including the 
empty ei) inside the background of ^/states. With the 
help of the fat-band representation, we fitted the rele- 
vant bands along the T-X direction (r-M in the usual 
two-formula unit cell) with parabolas shown as solid and 
dashed lines in Fig|5J Disregarding the possible warp- 
ing of the cylindrical Fermi surfaces, the DOS deriving 
from the above mentioned bands was estimated from the 



calculated effective masses by the free two-dimensional 
electron gas model N(E) — m* /irti 2 . This approxima- 
tion is even more justified by the fact that in the Eliash- 
berg analysis reported below only ratios between DOSs 
(which are much more accurate) enter in the calculations. 
Then, given this assumption, only the band curvature is 
relevant and therefore energy shifts with respect to LDA 
are not important in the model. The only important dif- 
ference relies on the fact that the number of parabolic 
bands taken into consideration is five or four. In this 
regard, the position of Ep is critical since band ei is 
just above Ep in the undoped compound. The super- 
conducting phase is obtained by 15% F doping i.e. in 
SmFcAsO0.85F0.15- F substitution adds electrons to the 
system and in the rigid band approach such effect can 
be coped by the rigid shift of Ep by 0.15 electrons per 
formula unit. Given the DOS N(E) of the system, the 
shift can be simply estimated as 0.15/Af (£p). Anyway, 
care should be taken in defining N(E), since the localized 
nature of Sm 4f orbitals makes it likely that they do not 
receive the additional doping charge from F. We there- 
fore filtered out this contribution from N(E) considering 
only the contribution of bands ei.2 and hi. 2. Following 
this approach we get a sort of upper bound for the Ep 
shift to be about 75 meV (90 meV) for SmFeAsO0.s5F0.15 
(SmFe0.75Ru0.25AsO0.85F0. 15)- The Fermi levels E F of F- 
doped and undoped SmFeAsO and SmFeo.75Ruo.25AsO 
are shown in Fig(5ta)-(c). We see that in the F-doped 
systems the band ei is always partially filled justifying 
our assumption to include its contribution in Eliashberg 
calculations. As for the effect of Ru substitution, by com- 
paring Fig[5Kb) and (c) we see that the effect is modest, 
only band 02 is a bit deeper and with lower effective mass 
as shown in FigJSJc). Calculated DOS and plasma fre- 
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x = 


m-eff 


DOS (st/Ha/Bohr 2 ) 


LJ P (meV) 


ei 


0.739 


0.2352 


714.04 


(-> 


1.508 


0.4800 


975.56 


hi 


0.908 


0.2890 


1391.91 


h 2 


2.057 


0.6547 


844.86 


x = 0.25 


ei 


0.668 


0.2127 


782.19 


f-2 


0.882 


0.2809 


1136.49 


fcl 


0.957 


0.3048 


1428.07 




1.736 


0.5526 


864.74 



TABLE II: Calculated effective masses, DOS and plasma fre- 
quencies for SmFei-^Ru^AsOo.ssFo.is at doping levels x = 
and x = 0.25. hi is doubly degenerate. 

quencies for each band and for the two doping levels con- 
sidered are reported in tabic HU 



Analysis of experimental results within Eliashberg 
theory 



• ab-initio 




FIG. 9: Temperature-dependent part of the superfluid den- 
sity, p s (T)/ p 3 (T = 0) (symbols) calculated by using the 
plasma frequencies obtained from first principles calculations 
(see table|lTJ in SmFei^Ru* AsO0.85F0.15 for x=0 (T c =52 K). 
The line is a fit with the formula tanh(l.76fc 1 / (T*/T - 1)), 
giving k — 0.6192. Inset: temperature-dependent part of the 
superfluid density (and therefore of the representative boson 
frequency [4(j) used in the Eliashberg calculations at different 
critical temperatures (see text for details). 



Based on the results of band structure calculations de- 
scribed so far, it is possible to propose an explanation 
of the experimental data by means of the simplest model 
that allows describing the essential physics of the materi- 
als under study. We used the three-band, s± Eliashberg 
theory [3l[ taking into account the feedback effect (ioj . 
Within this model we have two hole bands (from now 
on labeled as 1 and 2) and one equivalent electron band 
(labeled as 3). The free parameters are iV,(0), A31, A32, 
Q(T) and V. A^(0) is the DOS at Fermi level, calculated 
above for x = and x = 0.25 and obtained at all the 
other doping levels by linear interpolation as a function of 
the experimental T c \ A31 and A32 are the electron-boson 
coupling constants between band 3 and band 1 or 2, re- 
spectively; 51 (T) is the representative boson energy that 
we take as ft(T) = n tanh{1.76k^(T* /T - 1)) where 



Oo = 2T c /5 [4l|; the electron-boson spectral function 
has a Lorentzian shape [3l[ with halfwidth r = fio/2 @; 
T* is the feedback critical temperature, determined by 
solving the Eliashberg equations in the imaginary-axis 
formulation. The electron-boson coupling matrix is 

/ A 3 ii>i 
A 32 ^i 

y A31 A32 

where v x = JV 3 (0)/iVi(0) and v 2 = N 3 (0)/N 2 {0). k was 
determined at x = in the following way: first, the 
superfluid density was calculated by using the plasma 
frequencies obtained from first principles and reported 
in table HH giving also the correct T c . Then, since 
the temperature-dependent part of the superfluid density 
corresponds to that of the representative boson frequency 



[■101 . this curve was fitted with tanh(1.76fci/ (T*/T - 1)), 
giving k = 0.6192. The values of the calculated super- 
fluid density and the relevant fit are shown in Fig. |H] as 
symbols and line, respectively. 

At this point, the only parameters that remain to be 
determined are two: A31, A32. At x = they are ob- 
tained by reproducing the two experimental gaps, but 
it turns out that they also reproduce the exact experi- 
mental critical temperature. Then we assume that the 
ratio A31/A32 is doping- independent and that A: is a lin- 
ear function of the experimental critical temperature: 
k(T c ) = k(T c ^ x= o)T c /T CtX= o. This assumption is related 
to the temperature dependence of the superfluid density 
and is reasonable to suppose that with increasing x (and 
therefore decreasing of T c ) it decreases 0. Now, the only 
free parameter left in the process of fitting the T c depen- 
dence of the gaps is A31 which is fixed by obtaining T c 
coincident with the experimental one. The procedure is 
self-consistent since A31 is varied until T*, introduced in 
the formula for f2(T), allows reproducing the experimen- 
tal T c . Disorder effects have been neglected as impuri- 
ties are dominant in the intraband channel and are thus 
not pair-breaking. Moreover, we also assumed that, as a 
first approximation, they are also absent in the interband 
channel since the two gaps are well distinct at all doping 
levels. 

Results are shown in Fig. Htb) as solid lines: the cal- 
culated gap values follow rather well the trend of the 
experimental ones at all temperatures. The same good 
agreement with the experiment can be seen also by look- 
ing at the 2Ai/fcsT c ratios as a function of T c shown in 
Fig. [7] which also reports many other results from the 
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literature. In this regard, it is also remarkable that it 
is possible, with a relatively simple model and a small 
number of free parameters, to reproduce the increase of 
the ratio with decreasing T c . Since the strength of the 
coupling increases with decreasing T c , the same does, as 
expected, the total electron-boson coupling constant, X to t 
which is about 3.2 at T c =52 K and goes up to almost 7.3 
when T c =10 K. Another interesting result that comes 
out from the theoretical analysis is the temperature de- 
pendence of the representative boson frequency (shown 
in the inset to Fig. |H] for different critical temperatures) 
which, as already stated above, is equivalent to that of 
the superfiuid density [4(|. We can notice that, as T c 
decreases, the superfiuid density also decreases as a func- 
tion of T assuming, below T c =30 K, a positive curvature 
at intermediate temperatures. Similar dependencies (at 
least in the low- and mid-T c range) have been obtained 
in penetration depth measurements in Co-doped Ba-122 
samples, as reported in ref. Q. In that case the effect 
is less pronounced probably because the results are re- 
ported only down to about 2T c /5 while in our case T c 
drops to T c /5 at the highest Ru doping. Moreover, it is 
also interesting to notice that in Co-doped Ba-122 the 
temperature dependence of the superfiuid density looks 
slightly more depressed in samples that belong to the re- 
gion of coexistence of superconductivity and magnetism. 
This fact leads us to speculate that the observed behavior 
of the superfiuid density in our samples might be consid- 
erably influenced by the onset of a short-range magnetic 
order which competes with superconductivity and that 
has been observed by /iSR and 75 As NQR measurements 
in the same samples Penetration depth measure- 

ments, as the ones in ref. Q, would help clarifying this 
point as well as the experimental determination of the 
temperature dependence of the representative boson fre- 
quency, as done in refs. (3. EmL 142^ . 



CONCLUSIONS 

The isoelectronic substitution of Fe with Ru in opti- 
mally F-doped Sm-1111 is a good way to explore a very 
wide range of critical temperatures within the same Fe- 
based compound, in principle without changing the total 
charge of the system [l5|. The considerable decrease of 
T c induced by Ru substitution has been ascribed to dis- 
order in the Fe sub-lattice II 51 a nd /or to the onset of 
a short-range magnetic order [19J. Here we have shown 
that, in a wide range of T c (from 52 K down to 13.5 K, 
corresponding to Ru contents ranging from to 50 %) 
the Sm(Fei_ x Ru x )AsOo.85Fo.i5 system retains its origi- 
nal multi-gap character, and also the symmetry of the 
gaps remains nodeless. The amplitudes of the two ex- 
perimentally detectable gaps, Ai and A 2 , decrease al- 
most linearly with T c , but they remain well distinct down 
to the lowest T c . This suggests that the substitution- 



induced disorder mainly enhances intraband scattering 
and does not significantly affect the interband one. The 
gap ratios 2Ai/fcsT c strongly increase for T c < 30K in 
a manner which suggests an unexpected increase of the 
electron-boson coupling when T c is depressed. Very in- 
terestingly, the trend of the gap ratios as a function of 
T c in this single system is superimposed to the analogous 
trend obtained by plotting the data of many Fe-based 
compounds of different families [24|. Needless to say, 
this seems to point towards a general, universal prop- 
erty of this class of superconductors. By using the values 
of the density of states and of the plasma frequencies 
calculated from first principles, we have shown that the 
increase in the gap ratios 2Ai/ksT c on decreasing T c 
can be reconciled with a spin-fluctuation-mcdiatcd pair- 
ing even though the characteristic spin-fluctuation energy 
has been observed to decrease linearly with T c . The key 
to solving this puzzle is the feedback effect, i.e. the ef- 
fect of the condensate on the mediating boson, which is 
of course only expected when the superconducting pair- 
ing between electrons is mediated by electronic excita- 
tions [4^]. An analysis carried out within an effective 
three-band Eliashberg model shows indeed that the ex- 
perimental dependence of the gaps (and of the gap ratios 
2Ai/ksT c ) on T c can be explained as being due, in par- 
ticular, to a change in the shape of the temperature de- 
pendence of the characteristic boson energy, fi (T) , with 
respect to the optimal- T c compound (with no Ru). In- 
deed, a suppression of fio in the mid-temperature range 
(which becomes more and more sensible on decreasing 
T c ) is required to obtain the correct critical temperature 
and the correct gap values. This finding is in very good 
qualitative agreement with the experimental observation 
of a depression of the superfiuid density in Co-doped Ba- 
122 with reduced T c Q. The fact that in the latter case 
this reduction is observed in underdoped samples that 
fall in the region of coexistence of magnetism and super- 
conductivity further suggests that, also in our samples, 
the depression of the condensate (that in turns gives rise 
to a depression in the boson energy) at finite temperature 
may be considerably influenced by the onset of a short- 
range magnetic order competing with superconductivity 
induced by Ru substitution, as recently observed by fiSK 
and 75 As NQR measurements in the same set of samples 

El. 
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